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Topological effects and particle-physics analogies beyond the massless Dirac-Weyl 
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Armchair and zigzag edge terminations in planar hexagonal and trigonal graphene nanorings are 
shown to underlie one-dimensional topological states associated with distinctive energy gaps and 
patterns (e.g., linear dispersion of the energy of an hexagonal ring with an armchair termination 
versus parabolic dispersion for a zigzag terminated one) in the bands of the tight-binding spectra 
as a function of the magnetic field. A relativistic Dirac-Kronig- Penney model analysis of the tight- 
binding Aharonov-Bohm behavior reveals that the graphene quasiparticle in an armchair hexagonal 
ring is a condensed-matter realization of an ultrarelativistic fermion with a position-dependent mass 
term, akin to the zero-energy fermionic solitons with fractional charge familiar from quantum field 
theory and from the theory of polyacetylene. The topological origins of the above behavior are 
highlighted by contrasting it with the case of a trigonal armchair ring, where we find that the 
quasiparticle excitations behave as familiar Dirac fermions with a constant mass. Furthermore, 
the spectra of a zigzag hexagonal ring correspond to the low-kinetic-energy nonrelativistic regime 
of a leptonlike massive fermion. A onedimensional relativistic Lagrangian formalism coupling a 
fermionic and a scalar bosonic field via a Yukawa interaction, in conjunction with the breaking of 
the Z2 reflect ional symmetry of the scalar field, is shown to unify the above dissimilar behaviors. 

PACS numbers: 73.22.Pr, 73.21.-b, ll.lO.-z, 73.23.Ra 



I. INTRODUCTION 

Since its inception, relativistic quantum mechanics 
has been associated mainl;^ with the fields of particle 
and high-energy physics.^ ^ Recently, however, a table- 
top version of relativistic quantum physics emerged, fol- 
lowing the experimental isolation of graphene, which is a 
single-layer, planar honeycomb lattice of carbon atoms. ^ 
Indeed, the two-dimensional (2D) quasiparticle excita- 
tions of neutral graphene near the Fermi level behave^ 
as massless neutrinolike fermions described by the cel- 
ebrated Dirac- Weyl^ (DW) equation. The scientific 
and technological potential for exploiting charge carriers 
and quasiparticles with relativistic behavior in tunable 
condensed-matter and atomic-physics systems is attract- 
ing much attention.^^^^ In this context, an important 
question, as yet only partly explored, remains whether 
quasi-onedimensional (ID) graphene systems support ex- 
clusively DW massless or constant-mass Dirac fermions, 
or they can induce relativistic-quantum-field (RQF) be- 
haviors that require consideration of position-dependent 
mass terms, reflecting generalized underlying bosonic 
scalar fields. 

In this paper, we show that planar graphene nanor- 
ings do indeed exhibit a rich variety of physics, ranging 
from sophisticated RQF regimes to more familiar cases of 
constant-mass fermions (both in the relativistic and non- 
relativistic regimes). The emergence of these physical 
regimes depends on the specific combination of topologi- 
cal factors associated with modifications of the graphene 
lattice, such as the type of edge termination (i.e., arm- 
chair or zigzag) and the shape (i.e., hexagonal or trigonal) 
of the graphene ring. 

To this end, we investigate the properties of the cor- 



responding tight-binding^n^n^ (TB) spectra and of the 
associated Aharonov-Bohm^^ (the AB effect, which is a 
hallmark topological effect in condensed matter systems) 
oscillations of the magnetization, as a function of the 
magnetic field B. We then analyze the spectra and AB 
characteristics with the help of a Dirac-Kronig-Penney 
(DKP) superlattice model, ^ ^^ in the spirit of the vir- 
tual "band-structure" model for the nonrelativistic AB 
effect. ^^ We find that the relativistic behavior in arm- 
chair rings requires a profound modification of the ID 
Dirac equation^ through the introduction of a position- 
dependent mass term, in analogy with the fractional- 
charge, zero-energy topological modes in quantum field 
theory and in the theory of tran^-polyacetylene.^^^^ In 
contrast, the zigzag-ring spectra may correspond to the 
low-energy nonrelativistic regime of a leptonlike^^ mas- 
sive particle, heavier than the electron. 

Planar graphene rings (and the associated AB spectra) 
have been recently investigated by a number of groups us- 
ing tight-binding^^ 24 28 j^ethods (for polygonal shapes 
with armc hair or z igzag terminations), as well as con- 
tinuum DW *^^ * ^^ ^ equations supplemented with infinite- 
mass boundary conditions (for idealized circular shapes) . 
These earlier studies did not address the question of pos- 
sible analogies to ID quantum field theoretical models 
and particle physics. We stress that a prerequisite to 
raising and answering this question is the introduction 
by us of the virtual DKP superlattice model for the AB 
effect. For a review on recent experimental studies of the 
Aharonov-Bohm effect in graphene nanorings, see Ref. 
31, 

In addition to planar graphene rings, graphene 
nanoribbons (GNRs) are another class of related quasi- 
ID systems. GNRs have attracted substantially more 



attention than graphene rings and their study gave rise 
to a vast body of theoretical^^ ^^ and experimental^^ ht- 
erature. For graphene nanoribbons, it was found that 
a gap Ao may open at the Fermi energy, leading to an 
apparent analogy with the constant-mass Dirac fermion 
[see Ref. 32(c)]. As elaborated below, for an armchair 
nanoribbon (aGNR), this gap arises from the topology 
of the armchair edge which mimicks the dimerized do- 
mains (i.e., formation of Kekule unequal carbon bonds) 
in tran^-polyacetylene.^^^^^ Then in analogy with the 
scalar Z2 kink-soliton associated with the Peierls tran- 
sition in ^ran5-polyacetylene (or equivalently with the 
Z2 kink-soliton used in the Jackiw-Rebbi fermionic RQF 
modeP^) the qualitative features of the (fermionic) AB 
spectra of armchair graphene rings can be understood as 
resulting from an alternation (or lack of it) of two degen- 
erate dimerized domains associated with the arms of the 
graphene ring. We stress that the effective dimerization 
in the armchair GNRs and armchair graphene nanor- 
ings has a topological origin imposed by the presence 
of the armchair edges, while the dimerizati on in trans- 
polyacetylene is due to the Peierls instability.^^^'^ These 
two different underlying processes, however, lead to sim- 
ilar results that are characterized by the breaking of the 



ID Z2 refiectional symmetry (see in particular Sees. IV A 
and W\ below) . 

Our findings of quasiparticles in graphene with gen- 
eral position-dependent, and/or constant (rest), masses 
(unlike the massless neutrinolike quasiparticle in 2D 
graphene) is particularly interesting in light of increas- 
ing current interest in mass acquisition mechanisms, e.g., 
the Higgs mechanism in elementary particles^^ and 
condensed-matter physics .E^ESl 

The predicted unprecedented emergent unfolding of 
fundamentally distinct physical regimes, namely com- 
plex quantum-field theoretical ones versus nonrelativis- 
tic constant-mass ones, depending solely on the mate- 
rials' shape and edge termination is to date unique to 
graphene. It will be of great interest to test signa- 
tures of such regime-crossover experimentally for specifi- 
cally prepared graphene systems with atomic precision,'^ 
as well as to explore possible occurrence of such 
topological-in-origin physical behavior in other designer- 
Dirac-fermion artificial systems^ ^^- or nanopatterned ar- 
tificial graphene.'^ 

The plan of the paper is as follows. 

Sec. |TI| describes the Aharonov-Bohm tight-binding 
spectra for three characteristic planar graphene nanor- 
ings, i.e., an armchair hexagonal ring, an armchair trig- 
onal ring, and a zigzag hexagonal ring. 

Sec. |III| introduces the theoretical aspects of a rel- 
ativistic ID Dirac-Kronig-Penney model, based on the 
generalized Dirac equation. The DKP model describes 
the virtual superlattice associated with the Aharonov- 
Bohm effect. 



Sec. |IV| presents the DKP interpretation of the tight- 
binding spectra calculated in Sec. [llj The quasiparticle 
excitations in graphene nanorings are shown to exhibit 



behavior associated with quantum field theoretical mod- 
els for elementary particles beyond the massless Dirac- 
Weyl fermion. 

Sec. [V| discusses the full relativistic quantum field La- 
grangian formalism that underlies the DKP interpreta- 



tion elaborated in Sec. IV The Lagrangian formalism 
shows that the physics of quasiparticle excitations in pla- 
nar graphene nanorings relates to mass acquisition and 
formation of fermionic solitons. 



Finally, Sec. VI presents our conclusions. 



II. TIGHT-BINDING CALCULATIONS 

In this section, we will describe the TB spectra and 
corresponding AB magnetizations (as a function of the 
magnetic flux) for three characteristic cases of planar 
graphene rings, and specifically for a hexagonal arm- 
chair ring, a trigonal armchair ring, and a hexagonal 
zigzag ring (all three of similar dimensions). We note 
that the arms of the armchair rings studied here corre- 
spond to the class of p erfect armchair nanoribbons re- 
ferred to as metallicp^^^ 133jt35] ^j^gy exhibit a vanishing 
energy gap, Aq = 0, between the valence and conduction 
bands in tight-binding and continuum DW calculations. 
However, any perturbation (including the incorporation 
into a nanoring structure) may resullP^^^ in the open- 
ing of a gap Aq > 0. The metallic aGNRs have a width 
corresponding to Mw = 3/ — 1, / = 1,2,3,... carbon 
atoms. Counting along a zigzag line in the middle of the 
arm (away from the corners), the hexagonal and trigonal 
armchair rings studied in this papei^^ have A/W = 14. 

To determine the single-particle spectrum [the energy 
levels Si{B)] in the tight-binding calculations for the 
graphene nanorings, we use the hamiltonian 



Hr 



TB 



<id> 



I' J ^i ^j 



h.c. 



(1) 



with <> indicating summation over the nearest-neighbor 
sites z, j. The hopping matrix element 



Lij — Lij exp 



K^cJ^. 



ds ■ A(r) 



(2) 



where r^ and Vj are the positions of the carbon atoms i 
and j, respectively, and A is the vector potential associ- 
ated with the applied constant magnetic field B applied 
perpendicular to the plane of the nanoring. 

The AB magnetization of the graphene ring is given 
by 



M(^) = -S 



dEtot 
d^ ' 



where the total energy 
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FIG. 1. (Color online) (a) Part of the hexagonal graphene ring (2718 carbon atoms) with armchair edges. 1 and denote 
the short and long carbon dimers, respectively. Lengths in units of the graphene lattice constant ao — 0.246 nm. (b) TB 
single-particle spectrum as a function of the magnetic flux (magnetic field). Energies in units of the TB hopping-parameter 
t — 2.7 eV. The dashed black line denotes the Fermi level for A/" = 14 electrons. The arrows highlight the band gaps, (c) 
Magnification of the TB spectrum around the ^i band gap. ^i ~ 80 K, and thus it is easily detectable experimentally. 
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FIG. 2. (Color online) (a) TB total energies (sum over single- 
particle energies including spin) for A^ = 14 quasiparticles. 
(b) Correposnding TB magnetization (in units of the Bohr 
magneton). Energies in units of the TB hopping-parameter 

t = 2.7 eV. 



is given by the sum over all occupied single-particle en- 
ergies; the index a runs over spins. ^ = BS is the mag- 
netic flux through the area S of the graphene ring and 
$0 = hc/e is the flux quantum. 

The diagonalization of the TB hamiltonian [Eq. (IT])] 
is implemented with the use of the sparse-matrix solver 
ARPACK."^^ In calculating ^tot [see Eq. (|4|], only 
the single-particle TB energies with £i{B) > are 
considered. ^^ 

The TB results exhibit significant differences between 



the three cases studied here. These differences fall into 
two categories, namely, (A) same edge termination but 
different shape and (B) same shape but different edge 
termination. 



A. Hexagonal versus trigonal armchair rings 

The shape of the hexagonal graphene ring with arm- 
chair edge terminations considered here, as well as the 
corresponding TB results regarding the single-particle 
spectrum, the total energy, and the magnetization are 
displayed in Fig. fl] and Fig.^ as a function of ^/^o- The 
shape of the trigonal graphene ring with armchair edge 
terminations considered here, as well as the correspond- 
ing TB results regarding the single-particle spectrum are 
displayed in Fig. |3] 

Both TB energy spectra in Fig. [lib) (armchair 
hexagon) and Fig. p[c) (armchair triangle) are organized 
in braid bands separated by energy gaps. They exhibit, 
however, two main differences. The first concerns the 
composition of the braid bands, with the hexagonal ring 
exhibiting six-membered bands while the trigonal ring 
having three- membered bands. The sixfold and threefold 
groupings are a reflection of the Zq and Zs point-group 
symmetry of the hexagonal and trigonal shapes, respec- 
tively; these symmetries are fully taken into account by 
the DKP modeling in Sec. |ni| 

The second important difference between the TB en- 
ergy spectra in Fig. [ijb) and Fig. [sj^c) concerns the num- 
ber and nature of energy gaps. Specifically, two regular 
superlattice gaps Si and 82 are present in both cases; 
note the similar energy scale [a magnification of the re- 
gion around the ^1 gap is shown in Fig. [ijc)]. However, 
while a mass gap Aq is well developed for the trigonal 
ring [Fig. IsFc)], no corresponding Aq gap is present in 
the spectrum of the hexagonal ring, where the £ = 
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FIG. 3. (Color online) (a) Part of the trigonal graphene ring 
(2142 carbon atoms) with armchair edges, (b) Magnification 
of the corner of the trigonal ring shown in (a). 1 and denote 
the short and long carbon dimers, respectively. Lengths in 
units of the graphene lattice constant ao = 0.246 nm. (c) 
TB single-particle spectrum as a function of the magnetic 
flux (magnetic field). Energies in units of the TB hopping- 
parameter t — 2.7 eV. The two arrows denoted ^i and 82 
highlight band gaps. The arrow denoted by Ao indicates the 
opening of a gap at the Fermi level [e — 0) associated with 
generation of a rest mass M. Note that the Ao gap is absent 
in the TB spectrum of the hexagonal graphene ring in Fig. 



horizontal axis dissects (splits in half) the correspon ding 
sixfold braid band [Fig.mb)]. As we will show in Sec. |IV[ 
the gap Aq, for the case of the armchair triangle, is con- 
sistent with the physics of a massive (but still relativistic) 
Dirac fermion, while the dissecting of the £ = sixfold 
band, in the case of the armchair hexagon, is consistent 
with the formation of a fermionic soliton^^ built on a 
scalar Z2 kink soliton (precisely, a train of six fermionic 
solitons attached to successive Z2 kink/antikink solitons; 
see also Sec. W\. 



B. Armchair versus zigzag hexagonal rings 

As aforementioned, the TB results exhibit also signifi- 
cant differences between the cases of hexagonal rings with 
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FIG. 4. (Color online) (a) Part of the hexagonal graphene ring 
(2688 carbon atoms) with zigzag edges. 1 and denote the 
short and long carbon dimers, respectively. Lengths in units 
of the graphene lattice constant ao = 0.246 nm. (b) TB single- 
particle spectrum as a function of the magnetic fiux (magnetic 
field). Energies in units of the TB hopping-parameter t — 2.7 
eV. The dashed black line denotes the Fermi level for N — 20 
quasiparticles. The arrows highlight the band gaps. 



armchair and zigzag terminations. One such difference 
concerns the 5-dependence of the single-particle energies 
£(<!>), which is piecewise linear for the armchair case [Fig. 
l^b)], but piecewise parabolic for the zigzag case [Fig. 
4^b)]; this maintains also in the total energies £^tot(^) 
[Fig.j2Va) and Fig.lsja)]. For the AB magnetizations [see 
Eq. (fsf], this results in a characteristically different pro- 
file for the AB oscillations of M(^): step-staggeredlike 
in the armchair case [Fig. [2|b)] and sawtoothlike [Fig. 
|5|b)] in the zigzag case. The paraboliclike 5-dependence 
in the zigzag edge case^ is reminiscent of the spectra of 
a nonrelativistic ideal metal ring.^'^ In contrast (see be- 
low), the linear 5-dependence (in conjunction with the 
other features of the hexagonal armchair spectrum) can 
be associated with the fully relativistic regime of Dirac 
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FIG. 5. (Color online) (a) TB total energies (sum over single- 
particle energies including spin) for N = 20 quasiparticles. 
(b) Correposnding TB magnetization (in units of the Bohr 
magneton). Energies in units of the TB hopping-parameter 

t = 2.7 eV. 



fermions with position-dependent masses P^^ 

Both the armchair and zigzag single-particle spectra 
for hexagonal are organized in six-member braid bands 
separated by energy gaps. This sixfold grouping is a re- 
flection of the Zq point-group symmetry of the hexagonal 
rings. The energy gaps in the zigzag case are compara- 
ble to the width of the braid bands, and both the gaps 
and the widths of the bands increase with higher energy 
[see Fig. pFb)]; this is consistent with a nonrelativistic 
Kronig-Penney model. ""^^ In Fig. p[b), there are three en- 
ergy gaps labeled as ^o, ^i, and 62- Unlike the relativistic 
regime, in the nonrelativistic limit a gap around e = is 
unrelated to the particle mass, and for this reason we use 
the symbol 5o (with a lower-case S) instead of Aq as was 
the case in Fig. IsFc). We note that in the nonrelativistic 
regime the effective mass of the quasiparticle excitations 
is proportional to the inverse of the second derivative of 
the (approximately) parabolic spectra; see, e.g., Ref. [48l 



III. DIRAC-KRONIG-PENNEY 
SUPERLATTICE 

As was shown for a semiconductor ring using a non- 
relativistic superlattice approach, the AB single-particle 
spectrum exhibits energy gaps demarcating Bloch bands 
when a scatterer is placed on the ring.^^ In this con- 
text, the energy gaps that appear in Fig.jljb), Fig. [sj^c), 
and Fig. |4jb) indicate that the AB effect in polygonal 
graphene rings should be analyzed and modeled with the 
help of ID Kronig-Penney-type superlattices, with the 



corners of the polygons providing a generalized analog 
to the "scatterers" of Ref. [T8| Specifically, we consider 
a ID relativistic Dirac-Kronig-Penney model with unit 
cells built out of square potential barriers (developed in 
Ref. |2]in the context of the physics of quarks). 

The DKP model considered here is based on the ID 
generalized Dirac equation, which has the form: 



[E-Vix)]!"^ ^ihvra 



dx 



/3(^(x)^ = 0, (5) 



with vf being the Fermi velocity of graphene, which re- 
places the speed of light c; vp/c ^ 1/300. V{x) is an 
electrostatic potential and (j){x) is a bosonic position- 
dependent scalar field. We note that Ref. 2 uses (j){x) = 
Aivp + Vs{x)^ with A4 denoting the rest mass of a Dirac 
fermion (including the massless case) and the term Vg (x) 
being referred to as the Lorentz scalar potential. Omit- 
ting the last term on the left of Eq. ([5| reduces this equa- 
tion to the massless Dirac- Weyl^ one that underlies the 
majority of studies in planar graphene. 

The fermion field ^ is a twodimensional vector 



^ 






(6) 



where the subscripts u and / stand for the upper and 
lower component, respectively. The 2x2 Dirac matrices 
a and f3 can be^^ any two of the three Pauli matrices 



c^i 



1 

1 



(^2 







C^3 







(7) 



For example in the Dirac representation, one has a^ = 
(Ji and f3^ = as. /is the 2x2 identity matrix. We 
stress that the energy spectra of the DKP model are in- 
dependent of the specific representation used for a and 
p. Below we will often use the notation m(x), instead of 
^(x), to stress the fact that (j){x) can be considered as a 
position-dependent mass term. 

In the DKP modeling of the TB results, we assign 
to the nth side (n = l,...,Ns) of the polygonal ring 
a number (J) of square potential steps (regions) denoted 

as {V^^\m\^^)^ i = 1, . . . , J. We note again that the 

(n) 

electrostatic potentials V-^ ^ enter the ID Dirac equation 
[Eq. ([5] through the energy term {E — F)/^, while the 

(n) 

mass terms m- replace the scalar potential in the term 
/3(j){x)'^; as a result these two potentials lead to different 
physical behavior in the relativistic regime. 

The building block of the DKP model is a 2x2 wave- 
function matrix ft formed by the components of two inde- 
pendent 2x1 spinor solutions of the onedimensional first- 
order Dirac equation, ft plays^ the role of the Wronskian 
matrix W^^ used in the second-order nonrelativistic KP 
model; it is defined as follows at a point x of the unit cell 
(here we use the Dirac representation): 



^Kix) 



piKx 



iKx 



-Ke 



-iKx 



(8) 
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FIG. 6. (Color online) Schematic representat ion of the (V- ,m 



entering in the DKP calculations: (a) in Sec. IV A (b) in Sec. IV B (c) in Sec. IV C Each side of the polygon is divided in 
three lengths Li, L2, and L3. Boxes outside a polygon (online yellow) represent values mo > 0. Boxes inside a polygon (online 



1, . . . , 3 and n = 1, . . . , A/'s square-step parameters 



red) represent values —mo < 0. The nonzero values for the V^ and ^3 in (c) are portrayed by thick lines (online brown) in 
the interior of the schematic hexagon. Zero values of parameters are not highlighted. 



where 



K' 



2. ,4 



{E - Vy - m^v 



A 



hvpK 



E-V 



,,2 • 



We note again that, unlike the case of the original Dirac 
equation/ m here is not a constant, but it may take dif- 
ferent values from one region to the next. The trans- 
fer matrix for a given region (extending between two 
matching points Xi and X2) is the product Mx(2:i, X2) = 
Q,k{x2)0,'^{xi)] this latter matrix depends only on the 
width X2 — xi of the region, and not separately on xi or 
X2' The relativistic M matrices defined here correspond 
to those considered^^ in the case of a nonrelativistic su- 
perlattice in Ref. 14, 

The transfer matrix corresponding to the nth side of 
the hexagon is the product 



where we have explicitly denoted the dependence of the 
r.h.s. on the energy E. The presence {rj = 1/2) or ab- 
(9) sence {rj = 0) of an additional flux-shift in the relativis- 
tic or nonrelativistic case, respectively, follows through 
a comparison of the patterns of AB oscillations of a 
Dirac/Schrodinger electron in the limiting case of an ideal 
metallic circular ringP^l^ 



tn = ]^ MK{Xi,Xi^l) 
i=l,J 



xi = 0, xj+i = L, (10) 



with L being the (common) length on the hexagon side. 
The transfer matrix associated with the complete unit 
cell (around the ring) is the product 



Ns 



n=l 



where Ng is the number of sides of the polygonal shape 
considered {Ng = 3 and Ns = 6 for a triangle and 
hexagon, respectively). 

Following Ref. 18, we consider the supperlattice gen- 
erated from the virtual periodic translation of the unit 
cell as a result of the application of a magnetic field B 
perpendicular to the ring. Then the AB energy spectra 
are given as solutions of the dispersion equation 



cos[2^($/$o + r/)]=Tr[T(^)]/2, 



IV. DKP INTERPRETATION OF 
TIGHT-BINDING CASE STUDIES 

In this section, we demonstrate that our DKP model- 
ing can capture the essential physics underlying the TB 
spectra in Fig. [ijb). Fig. [sj^c), and Fig. |4jb). In this re- 
spect, it also provides a framework for unifying the broad 
variety of behaviors of the TB spectra of graphene nanor- 
ings. A schematic representation of the parameter sets 
used in our DKP simulations is given in Fig. [6] 



A. Armchair hexagonal ring and 
position-dependent-mass relativistic regime 

(11) First we attempt a solution corresponding to the gen- 
eralized Dirac equation ([5| with (j){x) = 0. In Fig. [Wa), 
we display the DKP spectra of a massless excitation (i.e., 

m-^ = for all i and n = 1, . . . , 6). This massless DKP 
spectrum does not exhibit any gaps and it is strictly lin- 
ear and periodic (with period <l>o) as a function of ^; this 
correlates with the behavior of a free massless fermion, 
as in the case of 2D graphene. We note that this gapless 
spectrum remains unchanged even when we consider in 
addition electrostatic potential steps . V- > 0, a fact 
that is a reflection of Klein tunneling.^^^^ 

(12) However, the TB spectrum in Fig.mb) exhibits energy 
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FIG. 7. (Color online) Spectra from the DKP model (relativistic regime) corresponding to the schematic case (hexagon) in Fig. 
p[a); see text and Eqs. ( |13|15| ) for the full set of parameters employed, (a) mo = and any Vi. Note the absence of band gaps 
due to the Klein paradox, (bj mo — 0.01t/v%. Note the similarity of the spectrum with that of the armchair graphene ring in 
Fig. Il[b). (c) mo — O.SOt/t'l^. The horizontal lines result from suppression of the AB oscillations due to strong localization. 
The dashed line indicates the energy zero. 



gaps (denoted as 5i and 82^ and highlighted by arrows), 
which require consideration of potential barriers in the 
DKP modeling. In the spirit of earlier investigations of 
real-space superlattices in 2D graphenep^^^^we consider 



first a constant mass m 



(n) 



A1 > 0, and alternating 



^i =Vq> {) and V-^^^ = steps in consecutive regions 



(see Sec. |III[ ). However, calculations with this choice show 
an opening of an energy gap at £ = 0, a fact that con- 
flicts with Fig. [ijb); in addition, it does not preserve the 
particle-hole symmetry of the TB spectra. 

A crucial feature of the TB armchair spectra in Fig. 
[lib) is the presence of zero-energy states (at half-integer 
fluxes). In order to capture this feature, and in light of 
our failed choices (see above), we attempt next to use 
a non-vanishing position-dependent scalar field (j){x) [de- 
noted also as m{x)] in Eq. ([s]), and set the electrostatic 
potential V{x) = 0. Recalling certain key elements in the 
theory of tmn^-polyacetylene pertaining to zero-energy 
solitonic modes, ^^ ^^ we employ in our DKP transfer- 
matrix solution of Eq. ([5| a scalar potential (j){x) of the 
form m{x) = —m{—x). Consequently, we divide each 
side of the hexagon in three parts ( J = 3) of length 



(n) 



and assign values 



'(n) 



6, L 



(n) 



VI 



(n) 



V, 



(n) 



V, 



(n) 



0, 



and 



(n) (n) 



0, 



4") = (-l)"mo. 



(13) 



(14) 



(15) 



Note that the index n = 1,...,6 here is numbering 
the sides of the hexagon, and thus overall the position- 
dependent mass term in our model is antisymmetric 
around each corner of the haxagonal ring. A schematic 



H 



representation of the above parameters [Eq. ([T3|— Eq. 
([l5|] is given in Fig. |6ja). 

Fig. IrFb) displays the DKP spectrum calculated with 
the dispersion equation (IT2| using the above parame- 
ter set with mo = O.Olt/v^^ and a = 8ao, b = IbaoW^ 
One observes that, in addition to the piecewise linear 
5-dependence, the DKP spectrum faithfully reproduces 
the two other central features of the TB spectrum [Fig. 
TJb)]: (i) the zero-energy states at half-integer values of 
^/^o and (ii) the opening at higher energies of energy 
gaps demarcating emerging sixfold braid bands. 

The behavior of each arm of the armchair hexagonal 
ring as a domain similar to the dimerized domains of the 
tran^-polyacetylene has a deeper physical reason, which 
can be revealed if one considers each arm of the hexagon 
as a perturbed armchair graphene nanoribbon. Indeed 
analytic expressions for the energy dispersion of the aG- 
NRs have been recently derived [see Refs. 32 (^b) and [35]: 
they have the form 



E{k) = ±\tie 



ika 



t2e 



-ikb\ 



(16) 



with k being the wave vector along the direction of 
the edge, a = ao/{2^/3) and b = ao/V^, with gq = 
0.246 nm being the lattice constant of graphene. ti = 
-2t cos[p7r / {^/w + 1)] + ^^1, P = 1,2,...,A/W and ^2 = 
—t + (5^2, with Sti, St2 denoting the perturbation away 
from a perfect aGNR. The dispersion equation (16) is 
similar to the tight-binding one describing a onedimen- 
sional chain of carbon atoms with bonds (hopping ma- 
trix elements) of alternating strength ti and ^2 (Kekule 
structure). The spectrum E{k) exhibits a mass gap 
Ao = 1^1—^2! at k = 7r/{a -\- b). This behavior is 
analogous to that of the linear-chain lattice TB model 
for tran^-polyacetylene; see Eq. (2.1) in Ref. |2l] In 
particular, when ti 7^ ^2, Eq. (16) describes a single 



dimerized domain breaking the ID reflect ional symme- 



try; when ti = ^2 (metallic aGNR), it describes a sym- 
metric chain of atoms, which preserves the reflect ional 
symmetry. We note that the factor underlying the for- 
mation of dimerized domains in tran^-polyacetylene is 
the Peierls in stabil ity incorporated in the Su-Schrieffer- 
Heeger modelP^^ The corresponding factor in armchair 
graphene nanoribbons and rings is topological in nature, 
i.e., it is a reflection of the lattice distortions of graphene 
due to the edge termination and the shape. 

Further insight can be gained through the observa- 
tion that the armchair edge by itself reflects the car- 
bon dimerization. Indeed it exhibits shorter dimers 
(denoted by 1) alternating with longer ones (denoted 
by 0); see Fig. [ija). Thus each arm of the graphene 
ring corresponds to one of two equivalent domains 
. . . , 1, 0, 1, 0, . . . or . . . , 0, 1, 0, 1, Following this no- 
tation and going around a given corner, one gets symbol- 
ically . . . , 0, 1, 0, 1, 1, 0, 1, 0, . . ., i.e., each corner (denoted 
by an underline) acts as a domain wall separating two 
alternative domains. 

In an infinite polyacetylene chain and at the position 
of the domain wall, a strongly localized fer mionic soli- 
ton develops having a fractional charge ±1/2.^^^^ How- 
ever, the graphene nanoring is a finite system and overall 
it can carry only integer charges. ^^^^ Furthermore, the 
TB results, and their DKP analog in Fig. JT^b), indicate 
a moderate extent of localization at the corners, result- 
ing from a non- negligible tunneling between the corners. 
Strong solitonlike localization (exhibiting a 1/6 fractional 
charge at each corner) can be achieved for larger values 
of mo. Indeed a large mo can localize a massless fermion, 
as is known from the (theoretical) trapping methodology 
in 2D graphene referred to as the infinite-mass bound- 
ary condition (introduced in Ref. 61 in the context of 
trapping neutrinos), as well as from the localization of 
massless quarks discussed in Ref. ^ The DKP spec- 
trum for a large value mo = O.SOt/vj^ is displayed in 
Fig. ItFc). It is seen that now the single-particle energies 
falling within the gap —0.3t < E < 0.3t form horizontal 
straight lines because the corresponding AB oscillations 
have been suppressed due to vanishing of the tunneling 
between the fermionic solitons at the corners; because 
of the localization no magnetic flux is trapped by the 
wave function on the hexagonal ring. It is of interest to 
note that such a train configuration of fermionic solitons 
in an hexagonal ring may be referred to as a fractional 
Wigner crystallite. We note that besides the zero-energy 
fermionic soliton discussed in the context of the states 
of polyacetylene,^^ ^^ Fig. mc) indicates the emergence 
of two polaroniclike state^^ with energies ~ ±0.16t 
falling within the gap. 



B. 



Armchair trigonal ring and constant-mass 
relativistic regime 



m 

Fig. IsFa) displays the DKP spectra of massless excita- 
tions (i.e., m^i' = for all i and n = 1,2,3). As was 
the case with the hexagonal ring [Fig.l7[a)], the massless 
DKP spectrum in Fig. [sja) does not exhibit any energy 
gaps and it is strictly linear and periodic (with period 
^o) as a function of <I> (or equivalently the magnetic field 
B). We again note that this gapless spectrum remains 
unchanged even when we consider in addition electro- 
static potential steps, l^ > 0, a fact that is a reflection 
of Klein tunneling. ^^ ^^ 

However, the TB spectrum in Fig. IsFc) exhibits en- 
ergy gaps (denoted as Aq, ^i and (^2, and highlighted 
by arrows), which require inclusion of potential barriers 
in the DKP modeling. Following the analogy with the 
tran^-polyacetylene, it is apparent that the opening of 
the Aq gap indicates the absence of domain alternation. 
Namely, the armchair trigonal ring represents a realiza- 
tion of a single domain extending along the full length 
of the triangle. Thus the corners of the triangle act a 
scatterers instead of domain walls as in the case of the 



hexagonal ring (see Sec. IV A). Indeed, following the no- 



tation of 1 and introduced above for the short and long 
dimers [Fig. Ilia)], and going around the corner of the 
trigonal ring in Fig. p[a) [or Fig. pfb)], one gets the se- 
quence . . . , 0, 0, 0, 1, 1, 1, 0, 0, 0, . . ., which is in agreement 
with the presence of the same domain on both sides of 
the 7r/3 corner. Consequently, in the DKP modeling we 
keep the same parametrization as in Eqs. (13) and (14), 
but we replace Eq. ( 15 ) by 



m 



(n) 



m 



(n) 



0, 



m 



(n) 



mo = M, 



1,2,3, (17) 



The results of our DKP calculations [associated with 
Eq. (p|] for the armchair trigonal ring are shown in Fig. 



that is, the mass parameters are the same (not alternat- 
ing) on the three arms of the trigonal ring [see schematic 
representation in Fig. [6]^b)]; recall that M denotes the 
fermion rest mass [see discussion below Eq. ([5|]. 

Fig. Sjb) [see also a magnification in Fig. [SJ^c)] dis- 
plays the DKP spectrum calculated from the dispersion 
equation (12) u sing the above parametrization [see Eqs. 
([l3|, ([m)), and ([l7|] with M = 0.02t/vj,, and a = 19ao, 
b = lOao-^^ One sees that the DKP spectrum reproduces 
the two central features of the TB spectrum in Fig. [sj^c) : 
(i) the gap Aq around e = and (ii) the threefold braid 
bands which are separated by the gaps ^i and 82 at higher 
energies. 



Zigzag hexagonal ring and constant-mass 
nonrelativistic limit 



As aforementioned, the TB spectrum of the hexagonal 
graphene ring with zigzag edges [Fig. [4Fb)] shows trends 
associated with a large-mass fermion in the nonrelativis- 
tic regime: (i) almost-perfect parabolic 5-dependence 
and (ii) energy gaps ^o, ^1 and 62 that are as large as 
the width of the braid bands. We remind that, for a 
massless excitation (ultrarelativistic limit), the spectra 
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FIG. 8. (Color online) Spectra from the DKP model (relativistic regime) corresponding to the schematic case (triangle) in Fig. 
p[b); see text and Eq. (17) for th e fu ll set of parameters employed. Panels (a) and (b,c) correspond to different choices of mo 
[denoted also as A^ here; see Eq. (17)]. (a) mo = and any Vi. Note the absence of band gaps due to the Klein paradox. (b,c) 
mo = 0.02t/vp, and in (c) we display a magnification of the region around E = 0. Note the similarity of the spectrum with 
that of the armchair graphene ring in Fig. ISTc). The dashed line indicates the energy zero. 
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FIG. 9. (Color online) Spectra from the DKP model (nonrela- 
tivistic regime) corresponding to the schematic case (hexagon) 
in Fig. Imc); see text and Eqs. ( 18][2Q ) for the full set of pa- 
rameters employed. The energies do not contain the rest- 
mass contribution p^ that is, E{t) — E — Mv%^ where E is 
the value obtained from solution of the DKP model and M 
is determined by fitting E to the TB results [Fig.|4];b)]. The 
rest mass determined to yield the best fit shown here, is given 
hy M — 42.06t/vF. This mass is heavier than that of the 
electron {2.10t/vp), suggesting an analogy with leptons. 



are strictly linear and exhibit no gaps [see Figs. iTFa) and 
[SJ^a)]; consequently small or moderate mass terms (rel- 
ativistic regime) will result in smaller energy gaps com- 
pared to the width of the braid bands [see Fig. W[h) and 
Figs, ^b) and gc)]. 

The sharply different physics (nonrelativistic versus 
relativistic behavior) underlying the TB spectra of the 
zigzag and armchair rings originates from the different 
edge topology. Following the notation of 1 and intro- 
duced above for the short and long dimers [Fig. ma)], and 



B 



going around the 27r/3 corner of the zigzag ring in Fig. 
H[a), one further gets the sequence . . . , 0, 0, 0, 1, 0, 0, 0, . . . 
This, in addition to the anticipation of a nonrelativistic 
regime, suggests that the entire hexagonal zigzag ring 
constitutes a single domain with embedded (corner) im- 
purities (electrostatic potential scatterers). 

In light of the above, using a large rest mass M in 
Eq. ([9| [nonrelativistic limit, see Eqs. (4.27) and (4.28) 
in Ref. 2 , we were able to reproduce in the DKP model 
[see Fig. [9] the overall trends of the TB spectrum [Fig. 
mh)] for the zigzag ring. The DKP parameters used^ 
were as follows. Lengths: 



4") = 

Rest mass: 
M 



-in) 



1.5an 



m 



(n) 



m 



(n) 



m 



-in) 



in) 



28ao. 



(18) 



42.06t/4. (19) 



Potential barrier step at each corner: 



V = V, 



in) 



V, 



in) 



80 X 10"^t, with V^ ^ = 0. (20) 



in) 



A schematic representation of these parameters is given 
in Fig.[6];c). 

In the evaluation of the nonrelativistic^ limit of the 
DKP model, one often sets E = Mvp + E for the pos- 
itive energies [see Eq. (4.27) in Ref. 2 . The quantity 
calculated from the DKP model in this limit is E^ while 
in comparing with the TB results [Fig. p[b)] we plot in 
Fig. [9|£^ vs. ^. The value of M. is determined by finding 
the best fit to the TB results.!^ It is expected that fur- 
ther improvement in the agreement between the TB and 
DKP approaches can be achieved by employing smooth 
(rather than square-shaped) profiles for the potential bar- 
riers. The constant mass, A^, found by us here is twenty 
times larger than the rest mass of the electron, indicat- 
ing an analogy with electronlike leptons, ^^ rather than 
the electron itself. 
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We note that the large rest mass M. found in this sec- 
tion is unrelated to the energy gap So [see Fig. [4]^b)]; 
{M.Vp/5o ~ 10^). Indeed in the nonrelativistic regime, 
the effective mass is related to the second derivative of 
the band spectra.^^ Our findings concerning the zigzag 
graphene ring are in agreement with the analysis of Ref. 
|48^ regarding the TB spectra of (infinite) symmetric poly- 
acene, which is a single chain of fused benzene rings 
and can be considered as the thinnest possible zigzag 
graphene nanoribbon. 



V. THE UNDERLYING RELATIVISTIC 
QUANTUM FIELD LAGRANGIAN 

In the above (see Sec. |ll|), we identified characteris- 
tic patterns of tight-binding spectra in planar graphene 
nanorings and we described their dependence on the edge 
termination (armchair versus zigzag) and the ring shape 
(hexagonal versus trigonal). Subsequently, we presented 
a unified interpretation of t he T B results using a Dirac- 
Kronig-Penney model (Sec. |III[ ), built upon the gener- 
alized ID Dirac equation [Eq. ([5|]. A central finding 
of our relativistic DKP analysis were the close analogies 
(found in Sees. IV A and |IVB ) between the behavior of 
armchair hexagonal and trigonal graphene nanorings and 
the physics of tran^-polyacetylene.^^ ^^ 

A natural next step towards a deeper understanding 
of the connection of our findings to relativistic quan- 
tum field theory is the elucidation of the underlying La- 
grangian formalism. Motivated by the theory of trans- 
polyacetylenej^^' we write a total Lagrangian density 



£ = £/ + £ 



(21) 



which is the sum of (partial) Lagrangian densities for the 
fermion field ^ [Eq. (|6|] and the bosonic scalar field ^(x). 
We note that the scalar field (j){x) was denoted earlier 
also as m{x) and was referred to as a position-dependent 



mass; see Sec. Ill and Sec. IV A 



The stationary generalized Dirac equation [Eq. ([5| 
be derived from the fermionic Lagrangian density 



can 



Cf = -ih^"^ — ^ - ihvF'^'^a — ^ - 0^^/3^, (22) 
at ox 

where we have neglected contributions from the electro- 
static potentials (see discussion in Sec. IV A). 
In Eq. (22), the last term 



jCy = -^^^/3^, 



(23) 



which depends on both the fermion ^ and scalar (j) fields, 
has the form of a Yukawa coupling. Cy is the potential 
agent for rest-mass acquisition by the originally massless 
fermion [described by the first two terms in the right 
hand side of Eq. ([22|]. 

The Yukawa interaction is also used in the Standard 
ModeP^"^ to describe the coupling between the Higgs 
field and the massless quark and lepton fields (i.e., the 



fundamental fermion particles). Through spontaneous 
symmetry breaking^^^^ of the Higgs field [which is a 
complex SU{2) doublet of four real scalar fields ^], these 
fermions acquire a mass proportional to the vacuum ex- 
pectation value of the Higgs field.'^^l'^ 

The essential observation that we make here is that, 
although, due to the ID character of the graphene rings, 
the ID Lagrangian in Eq. (21) does not possess the full 



richness of the Lagrangian of the Higgs sector in the 
Standard Model, both share the central aspect of sym- 
metry breaking and mass acquisition by a fermion via a 
Yukawatype interaction. 

We turn next to the task of constructing the La- 
grangian part C(j) for the scalar field ^, which is of the 
general form^^ ^^ (in 1 + 1 dimensions, i.e., time plus one 
space dimension) 



Ci 



IM2 



v{<P), 



(24) 



and which preserves the refiectional Z2 symmetry, i.e., it 
is invariant under cf) ^ —cf). Note that we are interested 
in the adiabatic approximation, and thus we omit the 
time dependent terms in iZ^. 

The emergence of a c onstant mass for the armchair 
trigonal ring (Sec. IV B) could simply be accounted for 
by considering a constant value of (j) = (j)o = M in the 
fermion Lagrangian Cf] then one needs to pay no further 
consideration to the bosonic £0. However, a more gen- 
eral position-dependent field, (j){x) = m(x), was found to 
be essential in our DKP-model analysis of the armchair 
hexagonal ring (Sec. IV A). In this case, (j){x) alternates 
between two unequal values ±(/)o, with (j)o = mo [see Fig. 
|6Fa)]. This indicates breaking of the Z2 symmetry of 
the solutions to the equation of motion derived from the 
Lagrangian in Eq. (24). An expression for V{(p) which re- 
produces qualitatively the above behavior (including the 
trigonal ring case) is the socalled 0^, which corresponds 
to a quartic double well potential in ^, i.e.. 



where ^ and ( are parameters. 



2\2 



C^) 



(25) 



With the potential in Eq. (25), the bosonic sector has 



the field equation (see Ch. 2.3 in Ref. [741 and Ch. 1.1 in 
Ref.[75|) 



dx'^ 



-^((/.2-c'),^ = o. 



(26) 



Two solutions of Eq. ([26]) are (j){x) = ±00 = ±C5 these 
solutions break the symmetry since 0o 7^ 0- Using these 
solutions in the Dirac Eq. ([5|, one obtains the standard 
constant-mass Dirac equation^ for the fermionic field. 
This case corresponds to the behavior of the armchair 
trigonal ring (Sec. IV B ), as well as to that of the zigzag 
hexagonal ring (Sec. IV C). 



In addition, however, Eq. (26) has nonlinear solutions 
that interpolate between the locations (j)o and —(po of the 
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FIG. 10. (Color online) Dashed line: The scalar Z2 kink 
soliton [Eq. (|27[)]. Solid line: The particle density (unnor- 
malized) of the corresponding enslaved fermionic soliton [Eq. 
(29)]. The values ^ = 1 and C — ^ were used. The domain 
wall is located at x = 0. 



two minima of the 1^(0) potential. One of these nonlinear 
solutions is called the Z2 kink soliton and the other the 
Z2 antikink soliton. The kink soliton is given by 




0/e(x) = Ctanh 

and the antikink soliton has the form 

^k{x) = -(t)k{x). 



(27) 



(28) 



We note that ^/c(±oo) = ±C and ^/c(±oo) = ^=("7 while 
0^(0) = 4(0) = 0. 

Using these kink or antikink scalar fields, the corre- 
sponding generalized Dirac equation [Eq. (Isl)] possesses 
fermionic soliton^^ solutions of the form (here a = cr2, 

/? = CTl) 



^ s{x) oc 







(29) 



The importance of these fermionic solitons lies in the 
fact that they have strictly zero energies\^^ thus they 
fall into the particle/antiparticle (valence/conductance 
band) energy gap. Furthermore, they are localized at 
X = 0, which is the domain wall between x > {(J)q) and 
X < (— ^o)- The enslavement of the fermionic soliton 
^ s{x) by the scalar potential of the kink soliton (j)k{x) 
is evident from Eq. (29). This is also reflected by the 
localization of ^ s{x) on the domain wall x = (see Fig. 

m. 



It is apparent that such zero-mode fermionic solitons 
^ s{x) underlie qualitatively the behavior of the fermion 
excitations in the armchair hexagonal graphene nanoring 
(Sec. IV A| ), with the corners of the ring behaving as do- 
main walls and the m{x) stepwise function in the DKP 
model [(j){x) in Eq. ([s])] mimicking an alternation of Z2 
kink and antikink scalar solitons [Eqs. (27) and (28)]; for 
a quantum-field-theory description of a train of alternat- 
ing kinks and antikinks, see Ch. 1.7 in Ref. [75l 



VI. CONCLUSIONS 

The paper investigated the different behavior of the 
Aharonov-Bohm spectra and magnetic-field induced os- 
cillations for three characteristic cases of planar graphene 
nanorings, i.e., an hexagonal ring with armchair edge 
terminations, a trigonal ring with armchair edge termi- 
nations, and an hexagonal ring with zigzag edge termi- 
nations. The tight-binding results (Sec. [ll|) were an- 
alyzed with the help of a ID relativistic Dirac-Kronig- 
Penney model^ (Sec. Ill), which accounts for the vir- 
tual superlattice associated^^ with the applied magnetic 
field. This analysis revealed unexpected topological ef- 
fects and condensed-matter analogies with elementary 
particle physics. 

In particular, the behavior found by us for the armchair 
hexagonal ring (Sec. IV A) is reminiscent of the extreme 
relativistic regime describing zero-energy fermionic soli- 
tons with fractional charge in quantum field theory^ and 
in the theory of tran^-polyacetylene.^^^^This regime re- 
sults from a consideration of a modified (generalized) 
Dirac equation with a position-dependent mass term (or 
equivalently a position-dependent scalar bosonic field). 
In contrast, the quasiparticle excitations in the arm- 
chair trigonal ring (Sec jIVB ) behave as relativistic Dirac 
fermions having a constant mass. A unification of these 
two dissimilar behaviors was presented in Sec. |V| by in- 
troducing the underlying relativistic Lagrangian formal- 
ism for a fermionic and a scalar bosonic fields coupled 
via a Yukawa interaction. The Yukawa term in conjunc- 
tion with the breaking of the Z2 reflectional symmetry 
of the scalar field may result in two outcomes, i.e., for- 
mation of a fermionic soliton (armchair hexagonal ring) 
or mass generation (armchair trigonal ring). The pro- 
foundly differing behaviors found by us for the armchair 
hexagonal and trigonal rings (with both sharing similar 
spatial dimensions), are manifestations of the quantum 
topological nature of this behavior, as distinguished from 
"quantum size effects" which are length-scale dependent 
phenomena, originating from spatial confinement of the 
electrons (quasiparticles, in general) .^^^^^^ 

The behavior of the zigzag hexagonal ring resembles 
the low-kinetic-energy nonrelativistic regime of a lepton- 
like fermion having a rest mass larger than that of the 
electron (Sec. IV C). This behavior contrasts with the 
relativistic ones found for the aforementioned armchair 
rings, thus highlighting the compounded topological and 
edge -termination effects. 

These findings^^ ^^ highlight the potential of graphene 
nanosystems for providing a bridge between condensed- 
matter and particle physics, well beyond the paradigm 
of the massless neutrinolike fermion familiar from the 
2D graphene sheet. Furthermore beyond the realm 
of graphene proper, where atomically precise narrow 
nanoribbons have already been synthesized,'^ we antic- 
ipate that our theoretical predictions could be tested 
using an ever expanding class of designer-Dirac-fermion 
manmade systems, such as optical lattices comprising ul- 
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tracold atoms j^^^^ or "molecular"^^ and nanopatterned 
artificial graphene™ 
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